Accepted for publication in the Ap.J. Letters 

Preprint typeset using LAT^X style emulatcapj v. 04/03/99 



O 



INTERVENING O VI QUASAR ABSORPTION SYSTEMS AT LOW REDSHIFT: 
A SIGNIFICANT BARYON RESERVOIR 1 

Todd M. Tripp, 2 Blair D. Savage, 3 and Edward B. Jenkins 2 

Accepted for publication in the Ap.J. Letters 

ABSTRACT 

Far-UV echelle spectroscopy of the radio-quiet QSO H1821+643 (z cm = 0.297), obtained with the Space 
Telescope Imaging Spectrograph (STIS) at ~7 km s _1 resolution, reveals 4 definite O VI absorption line 
systems and one probable O VI absorber at 0.15 < z a b s < 0.27. The four definite O VI absorbers are 
located near galaxies and are highly displaced from the quasar in redshift; these are likely intervening 
systems unrelated to the background QSO. In the case of the strong O VI system at z a b s = 0.22497, 
multiple components are detected in Si III and O VI as well as H I Lyman series lines, and the differing 
component velocity centroids and b— values firmly establish that this is a multiphase absorption system. 
, A weak O VI absorber is detected at z a b s = 0.22637, i.e., offset by ~340 km s _1 from the z a b s = 0.22497 

■ system. Lya absorption is detected at z a b s = 0.22613, but no Lya absorption is significantly detected 
at 0.22637. Other weak O VI absorbers at z a b s = 0.24531 and 0.26659 and the probable O VI system 
at 0.21326 have widely diverse O VI/H I column density ratios with N(0 VI)/A(H I) ranging from 
< 0.14 ± 0.03 to 5.2 ± 1.2. The number density of O VI absorbers with rest equivalent width > 30 mA 

— , in the H1821+643 spectrum is remarkably high, dN/dz ~ 48, with a high (90%) confidence that it is 

J> ■ greater than 17. We conservatively estimate that the cosmological mass density of the O VI systems 

C — | ^b(0 VI) ;> 0.0008/1^5 . With an assumed metallicity of 1/10 solar and a conservative assumption that 

' the fraction of oxygen in the O VI ionization stage is 0.2, we obtain £lh(0 VI) J> Q.004h^. This is 

| comparable to the combined cosmological mass density of stars and cool gas in galaxies and X-ray 

■ emitting gas in galaxy clusters at low redshift. 

Subject headings: cosmology: observations — galaxies: halos — intergalactic medium — quasars: 

■ absorption lines — quasars: individual (H1821+643) 

Q H - 1. INTRODUCTION of the Space Telescope Imaging Spectrograph (STIS) on 

' ' m ,. j ui 4. j? t • vi r\ • j-- HST, which provides a resolution of —7 km s _1 (FWHM) 

1 he resonance line doublet oi Li- like O VI is a sensitive . '. , , ,, T ,i • , 

«~ • , p, , „. . . . , , , ., with broad wavelength coverage. In this paper we present 

j~j , probe ot hot colhsionally ionized or warm very low density . «2 a d §3 e es Its o oe obable a d fo defi 

c/3 . photoionized gas in the intergalactic medium and galaxy ., ° „ ,„ °. . ,, tt,™, n, n 

(Ti . v i ti, o,, T ,,,„,, m Vno-7 co j , u , , , mte O VI absorption line systems m the S 1 IS HI 82 1+643 

^ halos. The O VI AA1031. 92, 1037.62 doublet has been de- , T n, ,. ,, • ,. ,. r „ , . , 

, j . i ,. , r^r<r\ -j r spectrum. In q4 we discuss the implications ot the high 

tcctcd m absorption toward QSOs over a wide range ot *\ „ , , , , . , ... * 

i , -r, , C1 . m • on r,nnn\ rn. 1 j. j rate of occurance of O VI absorbers at low redshift. Ihe 

redshifts (see hi m Iripp & Savage 2000). I he lowest red- ,. , . r . , • , , , 

u-£i. r\ in t_ i_ i- i i • j. u direct information we obtain about the highly ionized state 

shift O VI absorbers are particularly interesting because .,. „ ,, r^.,„ n , 

5-^ ' ii, j va- c i • iv. r\ar\ • uj. r v ot the gas from the presence of O VI allows us to estimate 

2 . the redshifts of galaxies near the QSO sight lines can be , ■ , , r ,i , ,„ , , , 

, , P , . , . , ,° i the baryomc content of these systems. We conclude that 

measured, and the relationship between the O VI absorber , , , , , . , r n 

,. , , j-jtt.j.1 U VI systems are likely to harbor an important traction of 

properties and environment can be studied. Furthermore, , , ,, , , 

,,. i • , i v i , , ,. , r the baryons at the present epoch, 

cosmological simulations predict that a substantial frac- J r r 

tion of the baryons in the universe are in a shock-heated 

phase at 10 5 - 10 7 K at low z (e.g., Cen & Ostriker 1999; 2. OBSERVATIONS AND ABSORPTION LINE 
Dave et al. 1999), and preliminary results indicate that MEASUREMENTS 
low-z O VI systems may indeed be an important baryon H1821+643 was observed with STIS for 25466 seconds 
reservoir (Tripp & Savage 2000). In a previous paper, Qn 199g June 25 with the medium reso lution FUV echelle 
Savage, Tripp, & Lu (1998) studied an intervening O VI mode (E140M ) and the .2x0.06" slit. 4 This STIS mode 
absorber associated with two galaxies at z « 0.225 in the prov ides a resolution of R = A/AA « 46000 or FWHM 
spectrum of the radio-quiet QSO H1821+643 using a com- w ? km g -i (Kimble et al 1998 ) The data were re duced 
bmation of low resolution Hubble Space Telescope (HST) ag described by Tripp & Savage (20 00) including the scat- 
spectra with broad wavelength coverage and a high reso- tered light correction developed by the STIS Instrument 
lution HST spectrum with very limited wavelength cover- Definition Team. The spectrum extends from -1150 to 
age. We have re-observed this QSO with an echelle mode mQ A with four small gapg betw een orders at A > 1630 

Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astronomy. Inc., under NASA contract NAS 5-2555. The data were reduced with the STIS 
Team software, and the research was supported by NASA through grants GO-08165.01-97A and GO-08165.02-97A from STScI. 
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A. Throughout this paper wavelengths and redshifts are 
heliocentric. 

We first searched the spectrum for O VI absorbers by 
checking for lines with the velocity separation and rela- 
tive line strengths expected for the doublet. This iden- 
tified the four definite O VI systems. We then searched 
for O VI lines associated with known Lya absorbers, and 
this revealed the probable system (see below). A selected 
sample of the spectrum is shown in Figure 0. This por- 
tion of the spectrum shows the O VI doublet at z a b s = 
0.22497 as well as a much weaker O VI doublet at z a b s 
= 0.22637. Both of these O VI absorbers are discussed 
in §3.1. In addition to the O VI systems in Figure [j], 
the STIS echelle spectrum shows new O VI absorbers at 
z&hs = 0.24531 and 0.26659 which have small equivalent 
widths; these systems are briefly discussed in §3.2 along 
with a possible O VI system at z a b s = 0.21326. The z a b s 
= 0.21326 system is a strong H I Lya/Ly/3 absorber with 
a > 4<T line detected at the expected wavelength of O VI 
A1031.93. However, the corresponding O VI A1037.62 line 
is blended with Milky Way S II A1259.52 absorption due 
to two high velocity clouds, so we consider this a probable 
but not definite O VI detection. The component structure 
establishes that this blend is mostly due to Milky Way 
S II, but it is possible that an O VI A1037.62 line of the 
right strength is present as well. In principle this could be 
proved by comparing the S II A1259.52 line strengths to 
the S II AA1250.58, 1253.81 line strengths. However, this 
does not yield a clear result at the current S/N level due 
to ambiguity of the continuum placement near 1259 A. 

Rcstframe equivalent widths (W r ) of absorption lines de- 
tected in the O VI systems, measured using the software of 
Scmbach & Savage (1992), are listed in Table |[ Note that 
the quoted errors in equivalent width include contributions 
from uncertainties in the height and curvature of the con- 
tinuum as well as a 2% uncertainty in the flux zero point. 
Integrated apparent column densities (Savage & Sembach 
1991) are also found in Table [ljwith error bars including 
contributions from continuum and zero point uncertain- 
ties. To measure line widths, we used the Voigt profile 
fitting software of Fitzpatrick & Spitzer (1997) with the 
line spread functions from the Cycle 9 STIS Handbook. 

3. ABSORBER PROPERTIES 

Four of the five absorption systems in Table [j] are within 
a projected distance of 1 Mpc or less of at least one 
galaxy with | Av | =| c(z ga i - z abs )/(l + z mC an) | < 300 
km s , and some of them are close to multiple galax- 
ies (see Table 1 in Tripp et al. 1998). These absorbers 
are also displaced from the QSO redshift by 7100 km 
s- 1 (z abs = 0.26659) to 17100 km s" 1 (z abs = 0.22497). 
Finally, the O VI profiles are relatively narrow. There- 
fore these are probably intervening systems that trace the 
large-scale gaseous environment in galaxy envelopes and 
the IGM rather than "intrinsic" absorbers (Hamann & 
Ferland 1999). 

3.1. VI Absorbers at z = 0.22497 and 0.22637 

Since the O VI doublets at z abs = 0.22497 and 0.22637 
shown in Figure |l| are separated by only ~340 km s _1 , 
they are probably related and we discuss them together. 

5 In this paper, the cosmological parameters are set to Hq = 75/175 



Two emission line galaxies are known at heliocentric red- 
shifts of 0.22560 and 0.22650 at projected distances of 105 
and 388 kpc from the sight line 5 (Tripp et al. 1998). 
In addition to the O VI doublet, the STIS spectrum shows 
strong absorption lines due to H I Lya , Ly/3, Ly7, Si III 
A1206.5, and possibly C III A977.02 at z abs = 0.22497; the 
absorption profiles of most of these species are plotted on 
a velocity scale in Figure 0. The lines of N V and Si IV 
are not detected at greater than 3<7 significance, and upper 
limits on their equivalent widths and column densities are 
listed in Table [l] along with upper limits on C II and Si II. 

Figure || provides several indications that these systems 
are multiphase absorbers. Several components are readily 
apparent in most of the z a b s = 0.22497 profiles including 
the O VI lines (see also Figure Ej). Fitting of the Si III 
profile yields b = 7.7j\j '? an d l-6_o' g km s _1 for the two 
well-detected components at v = -3 and +25 km s _1 , re- 
spectively. However, the component velocity centroids and 
6-values are not compatible with a homogeneous mixture 
of O VI and Si III. For example, the Si III profile shows a 
prominent narrow component at v « 25 km s _1 , and there 
is no obviously corresponding component in the O VI pro- 
files. While thermal Doppler broadening can make the 
O VI profiles broader than those of Si III, at most the in- 
crease will be a factor of y / 28/16, and this is inadequate 
to produce the breadth of the observed O VI lines. Thus 
we are compelled to consider a mixture of phases, some 
of which show up in Si III, while others are prominent in 
O VI. 

In the case of the O VI at z a bs — 0.22637, which is also 
visible in Figure g the multiphase nature is suggested by 
an offset of 60 km s _1 between the H I Lya and O VI 
velocity centroids. Also we note that no H I absorption is 
significantly detected at the velocity of the O VI suggesting 
that the hydrogen is thoroughly ionized in the O VI gas. 
This O VI absorber may be analogous to the highly ion- 
ized high velocity clouds seen near the Milky Way which 
show strong high ion absorption with very weak or absent 
low ion absorption (Sembach et al. 1999). 

3.2. Other Weak O VI Systems 

The two new O VI systems at z a b s = 0.24531 and z abs 
= 0.26659 are plotted in Figure ||. A striking feature of 
these weak O VI absorbers (and the candidate O VI at 
Zabs = 0.21326) is that while their O VI column densities 
are comparable, the strengths of their corresponding H I 
absorption lines are significantly different (see Table [l] 
and Figure g). For example, N(0 VI)/JV(H I) = 5.2±1.2 
in the z abs = 0.24531 system while N(0 VI)/AT(H I) = 
1.2±0.2 in the z a b s = 0.26659 absorber. The contrast 
is even more dramatic with the z a b s = 0.21326 absorber, 
which has N(0 Vl)/iV(H I) < 0.14 ± 0.03. For reference, 
in collisional ionization equilibrium (Sutherland & Dopita 
1993), gas with solar metallicity at the peak O VI ion- 
ization temperature should have N(0 Vl)/iV(H I) ~ 100. 
The large variablility of the observed O VI/H I ratio could 
indicate that the metallicity of the O VI absorbers varies 
substantially, or this could be due to differences in the 
physical conditions and ionization of the gas. If, for ex- 
ample, these are multiphase absorbers with the H I lines 
arising in a cool phase which is embedded in a hot phase 

km s -1 Mpc -1 and q = 0.0. 
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(Tripp & Savage 2000), and (2) one or two additional in- 
tervening O VI absorbers are evident in the H1821+643 
spectrum which did not satisfy the selection criteria to be 
included in the sample. More observations are needed to 
build the sample of weak O VI lines at low z. 

For comparison, low to moderate redshift Mg II ab- 
sorbers with W r > 20 ml have dN/dz = 2.65 ± 0.15 
(Churchill et al. 1999; see also Tripp, Lu, & Savage 1997). 
The stronger O VI absorbers are less common; Buries & 
Tytler (1996) report dN/dz = 1.0 ± 0.6 for O VI systems 
with W r > 210 mA at < z abs > = 0.9. Evidently, the 
dN / dz of the weak O VI lines is substantially larger than 
dN/dz of other known classes of low z metal absorbers and 
is more comparable to that of low z weak Lya absorbers, 
which have dN/dz ~ 100 for W r > 50 mA (Tripp et al. 
1998; Penton et al. 2000). 

Following analogous calculations (e.g., Storrie-Lombardi 
et al. 1996; Buries & Tytler 1996), 8 the mean cosmological 
mass density in the O VI absorbers, in units of the current 
critical density p c , can be estimated using 



n b {o vi) 



p c cf(0 VI) 



O VI 



EjjVjOg vi) 

AX 



(1) 



which produces the O VI absorption (e.g., Mo & Miralda- 
Escude 1996), then the wide variations in the O VI/H I 
ratio could simply be due to interception of fewer cool 
phase clouds in one absorber compared to another. 

A full analysis of the range of physical conditions of 
these absorbers will be presented in a later paper. How- 
ever, it is interesting to note that the H I Lya profile 
of the z a bs = 0.26659 system is rather broad and rela- 
tively smooth (see the bottom panel of Figure |J) , which 
may indicate that this absorber is collisionally ionized and 
hot. However, fitting a single component to the z a bs = 
0.26659 Lya profile yields b = 44.6lg'g, which implies that 
T < 1.2 x 10 5 K. At this temperature the O VI ionization 
fraction is rather small in collisional ionization equilibrium 
(Sutherland & Dopita 1993), and an unreasonably high 
metallicity is required to produce the observed A^(0 VI) 
and N(H I) in the same gas. This may be another indica- 
tion that these are multiphase absorbers or that the gas is 
not in ionization equilibrium. 

4. NUMBER DENSITY AND COSMOLOGICAL MASS 
DENSITY 

The new STIS data in this paper provide an opportu- 
nity to evaluate the number density of low-z O VI ab- 
sorbers per unit redshift (dN/dz) and a lower limit on 
their cosmological mass density. If we neglect continuum 
placement uncertainty and other systematic error sources, 
the STIS E140M spectrum of H1821+643 is formally ade- 
quate for > 4cr detection of narrow lines with W r > 30 mA 
at A obs £ H88 A (z abs £ 0.151 for O VI A1031.93). How- 
ever, the continuum placement ambiguity substantially in- 
creases the uncertainty in Wr for weak lines. Moreover, 
broader resolved lines spread over more pixels have higher 
limiting equivalent widths (limiting W oc y/no. pixels), so 
broad weak lines may not be detected at the 4tr level. Con- 
sequently, the dN/ dz derived below should be treated as a 
lower limit. We require detection of both lines of the O VI 
doublet with W T > 30 mA, and we exclude one absorber 6 
within | Av | < 5000 km s _1 of z cm to avoid contamina- 
tion of the sample with intrinsic absorbers. This results 
in a sample of three O VI systems 7 over a redshift path of 
Az — 0.063 (after correction for a loss of Az — 0.061 for 
spectral regions in which either of the O VI lines is blocked 
by ISM or extragalactic lines from other redshifts). There- 
fore the most probable dN/dz ~ 48 for W r > 30 mA and 
0.15 < z a bs < 0.27, and we conservatively conclude that 
dN/dz > 17 at the 90% confidence level (following the 
Gehrels 1986 treatment for small sample statistics). This 
is a remarkably high number density. It is important to 
emphasize that the sample is extremely small and, since 
very little is known about weak O VI lines at low redshift, 
it remains possible that dN/dz is unusually high toward 
H1821+643 for some reason. However, there is supporting 
evidence that dN/dz is generally high: (1) a similar dN/dz 
is derived from STIS echelle spectroscopy of PG0953+415 

6 We exclude the associated O VI absorber at z a b s = 0.2967. This system is not listed in Table h but is discussed in detail in Savage et al. 
(1998) and Oegerle et al. (2000). 

7 The three systems include those at z a t, s = 0.22497, 0.24531, and 0.26659. We exclude the probable system at z a b s = 0.21326, and the z a ^ s 
= 0.22637 system falls below the equivalent width threshold. 

8 Note that while Buries & Tytler (1996) calculated the cosmological mass density of the oxygen ions in O VI absorbers (which is quite 
small), they did not apply an ionization and metallicity correction to estimate the total baryonic content of the O VI systems. Instead, they 
used this method to place a lower limit on the mean metallicity of the O VI systems. 

9 For a discussion of the impact of small number statistics on the f2j,(0 VI) estimates, see Tripp & Savage (2000). 



where fi is the mean atomic weight (taken to be 1.3), 
/(O VI) is a representative O VI ionization fraction, 
(0/H)o vi is the assumed mean oxygen abundance by 
number in the O VI absorbers, J^i (O VI) is the 
total O VI column density from the i absorbers, and 
AX is the absorption distance interval (Bahcall & Pee- 
bles 1969), corrected for blocked spectral regions. With 
the sample defined above, we have f2b(0 VI) = 8.0 x 
10~ 5 /(O VLr 1 10- [°/ H l /if/ where [O/H] = log (O/H) - 
log (O/H) . To set a conservative lower limit on ftb(0 VI), 
we assume [O/H] = -0.3 and /(O VI) = 0.2 (which is 
close to the maximum value in photo- or collisional ion- 
ization, see Tripp & Savage 2000), which yields f2b(0 VI) 
> 0.0008/175 1 . If we set the mean metallicity to a more 
realistic value such as [O/H] = —1, flb(0 VI) increases 
to > 0.004/iyg 1 . Similar lower limits on fi;,(0 VI) have 
been derived by Tripp & Savage (2000) using a slightly 
less sensitive sample based on STIS echelle spectroscopy of 
PG0953+415 and earlier Goddard High Resolution Spec- 
trograph observations of H1821+643. The lower limit as- 
suming (O/H) = 1/10 solar is comparable to the combined 
cosmological mass density of stars, cool neutral gas, and 
X-ray emitting cluster gas at low redshift, f^-l-Jln i 2icm + 
^h 2 +^x-ray ~ 0.006 (Fukugita, Hogan, & Peebles 1998). 
Though still uncertain due to the small sample, 9 small red- 
shift path probed, and uncertain (0/H)o vi) these prelim- 
inary lower limits on Slj, (O VI) suggest that O VI absorbers 
contain an important fraction of the baryons in the low 
redshift universe. 
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Fig. 1. — Portion of the STIS E140M spectrum of H1821+643 showing the strong O VI absorption lines at 2 abs = 0.22497 and the weaker 
O VI absorber at z a b s = 0.22637. The calibrated flux is plotted vs. observed heliocentric wavelength, and the solid line near zero is the la- 
flux uncertainty. The line at 1266.9 A is an unrelated C III line from the absorption system at z a t, s = 0.2967. In this figure, the spectrum has 
been binned 2 pixels — » 1 pixel for display purposes only (all measurements in the text were made using the unbinned full resolution data). 
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Table 1 

Equivalent Widths and Integrated Column Densities 



\ , a 


Species 


Ao b 




AT d 
JV a 






(A) 


(mA) 


(^1U cm j 






z abs = 0.21326 




1474.93 


H I 


1215.67 


471±14 


> 19.4° 


1244.47 


H I 


1025.72 


134±9 


24.7 ±2.0 


1252.07 


O VI f 


1031.93 


38±9 


3.55 ±0.81 






z abs = 0.22497 




1488.97 


H I 


1215.67 


852±22 


> 34.4 C 


1256.41 


H I 


1025.72 


503±14 


> 155 c 


1191.27 


H I 


972.54 


340±17 


> 266 c 


1264.09 


VI 


1031.93 


185±9 


19.9 ± 1.2 


1271.05 


O VI 


1037.62 


110±10 


21.0 ±2.0 




Si II 


1260.42 


< 47 g 


<0.3is 


1477.93 


Si III 


1206.50 


108±9 


0.73 ±0.07 




Si IV 


1393.76 


48±23 h 


<1.3 h 




C IP 


1036.34 


<40 g 


<3.6 


1196.79 


c m?j 


977.02 


319±16 


> 8.9 d 




N V 


1238.82 


<52e 


<2.5 g 






z abs = 0.22637 




1490.57 


H I 


1215.67 


168±15 


3.7 ±0.3 


1265.53 


VI 


1031.93 


25±5 


2.4 ±0.5 


1272.49 


O VI 


1037.62 


21±5 


3.7 ± 1.0 






z abs = 0.24531 




1513.86 


H I 


1215.67 


45±9 


1.0 ±0.2 


1285.07 


O VI 


1031.93 


55±6 


5.2 ±0.6 


1292.15 


O VI 


1037.62 


39±6 


7.0 ± 1.1 






z abs = 0.26659 




1539.75 


H I 


1215.67 


177±12 


4.2 ±0.3 


1307.04 


O VI 


1031.93 


55±8 


5.1 ±0.8 


1314.26 


O VI 


1037.62 


32±8 


5.8 ± 1.3 



a Observed vacuum Heliocentric wavelength of the line 
centroid. 

b Restframe vacuum wavelength from Morton (1991). Os- 
cillator strengths used for these measurements were also ob- 
tained from Morton (1991) except for the Si II A1260.42 
/—value, which is from Dufton et al. (1983). 

c Restframe equivalent width integrated across all compo- 
nents. 

d Intcgrated apparent column density 7V a = J Na,(v)dv. 

e Saturatcd absorption line. Lower limits are derived from 
integrated apparent column densities with pixels with flux 
< set to their 3<r upper limits. 

f This line is probably the stronger line of the O VI dou- 
blet. However, the line identification is less secure than the 
other O VI lines in this table because the weaker O VI line 
at 1037.62 A is blended with Milky Way S II A1259.52 ab- 
sorption. With higher signal-to-noisc, it should be possible 
to confirm or refute this identification. 

S 4(T upper limit. 

h A 2(7 feature is detected at the expected wavelength, 
but we do not consider this significance adequate to claim a 
reliable detection. Consequently, we set an upper limit on 
N a which is the measured column density of the marginal 
line + 2(7. 

'The somewhat stronger C II A1334.53 line falls in a gap 
between orders. 

J The strength and velocity extent of this line is surprising 
compared to the strength and velocity extent of the Si III 
A1206.5 line. However, it is difficult to find convincing al- 
ternative identifications of this strong feature. The line may 
be a blend. 
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Fig. 2. — Profiles of absorption lines detected in the O VI system at z a b s = 0.22497, plotted versus restframe velocity where v — at 2 a b s 
= 0.22497. The tick mark at v ~ 340 km s 1 indicates the velocity of the £ a bs = 0.22637 absorber in the z abs = 0.22497 restframe. The 
dotted lines show the continua adopted for absorption line measurements, and the grey lines show the flux zero levels. The data in this figure 
are shown at full resolution (no binning has been applied). Note that the strong lines at v < —200 and > 330 km s _1 in the H I Ly7 panel 
arc due to the ISM Si II 1190.42 and 1193.29 A lines, respectively. 
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Fig. 3. — Absorption profiles of lines detected in the O VI absorbers at ,z a b s = 0.24531 (upper three panels) and at z a b s = 0.26659 (lower 
three panels), plotted versus restframe velocity. 



